Tunability of the dielectric response of epitaxially strained SrTiO,3 from first principles 
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The effect of in-plane strain on the nonlinear dielectric properties of SrTiOa epitaxial thin films 
is calculated using density-functional theory within the local-density approximation. Motivated by 
recent experiments, the structure, zone-center phonons, and dielectric properties with and without 
an external electric field are evaluated for several misfit strains within ±3% of the calculated cubic 
lattice parameter. In these calculations, the in-plane lattice parameters are fixed, and all remaining 
structural parameters are permitted to relax. The presence of an external bias is treated approxi- 
mately by applying a force to each ion proportional to the electric field. After obtaining zero-field 
ground state structures for various strains, the zone-center phonon frequencies and Born effective 
charges are computed, yielding the zero-field dielectric response. The dielectric response at finite 
electric field bias is obtained by computing the field dependence of the structure and polarization 
using an approximate technique. The results are compared with recent experiments and a previous 
phenomenological theory. The tunability is found to be strongly dependent on the in-plane lattice 
parameter, showing markedly different behavior for tensile and compressive strains. Our results are 
expected to be of use for isolating the role of strain in the tunability of real ultrathin epitaxial films. 

PACS numbers: PACS: 77.55. +f, 77.84.Dy, 77.22. Ej 
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I. INTRODUCTION 

Microwave dielectric materials strongly tunable by the 
application of an electric bias field are increasingly impor- 
tant for a variety of applications in microwave electron- 
ics, including tunable capacitors and oscillators, phase 
shifters, and delay linesi 2 * 3 " The dependence of dielec- 
tric constant e on bias field £ is typically strongly non- 
linear, and a tunability parameter n can be defined as 
n = [e(0) — e(£)]/e(0), where £ is an operational bias field 
of interest. Because of the potential technological im- 
pact, considerable experimental effort has been directed 
toward the design, control, and optimization of highly 
tunable low-loss materials. 

For many applications, thin-film morphologies are re- 
quired, and in recent years the growth of high-quality ul- 
trathin perovskite films with unprecedented atomic-level 
control has become possible, using techniques such as 
molecular-beam epitaxy (MBE)4 and pulsed-laser depo- 
sition (PLD)i These efforts have achieved the growth of 
single-crystal films of nanometer-scale thickness with a 
minimum of defects. Even so, the dielectric properties of 
these films often differ quite substantially from bulk 6-7 . 
This difference is due to a number of factors, of which 
the strain in the film is in many cases among the most 
important. More specifically, the mechanical boundary 
conditions on a coherent epitaxial film require the in- 
plane lattice constant of the film to stretch (or contract) 
to match the lattice constant of the substrate. Both 
experimental and theoretical studies have found that 
even small epitaxial strains can appreciably influence the 
Curie temperature and the dielectric permittivity and 
tunability of BaTi0 3 , SrTi0 3 , and Bai-^Sr^TiOa (BST) 
thin films, 8 ! 9 . 10 ! 11 

Paraelectric BST at x near 0.5 is already known to 



possess a large dielectric response and high tunability at 
room temperature, associated with close proximity to the 
ferroelectric phase transition at about — 23 K (Ref . IT^) . 
However, significantly lower loss has been reported for 
single-crystal SrTi0 3 (STO) thin films, 13 possibly due to 
the absence of compositional disorder. The ground state 
of STO is nonpolar but nearly ferroelectric, and thus 
small applied stresses would be expected to have a signif- 
icant influence on the Curie temperature and associated 
susceptibilities. Indeed, recent experimental studie s 14 ! 15 
indicate that significant changes in Curie temperature 
and dielectric properties occure when STO is grown epi- 
taxially on substrates with different lattice constants. 

One of the best-studied perovskite materials, pure 
STO is the subject of considerable experimental and the- 
oretical literature. STO adopts the centrosymmetric cu- 
bic perovskite structure at room temperature, and un- 
dergoes a structural phase transition from the cubic to 
a tetragonal, nonpolar antiferrodistortive (AFD) phase 
when cooled below 105 This transition, however, 

is observed to have little effect on the dielectric prop- 
erties. Cooling to still lower temperatures results in a 
strong Curie- Weiss-like increase in the static dielectric 
response, suggestive of a phase transition at about 20 K. 
However, no transition actually occurs in that tempera- 
ture range; instead, the dielectric constant saturates to a 
value of ~ 2 x 10 4 at zero temperature. 17 ' 18 ! 19 ! 20 ! 21 ! 22 ! 23 
The absence of the ferroelectric (FE) transition can be 
explained by quantum fluctuations of the atoms about 
their centrosymmetric cubic positions (i.e., the forma- 
tion of a "quantum paraelectric" state) i 19 i 20 i 21 i 22 i 24 i 25 
The proximity to a ferroelectric transition is also ev- 
ident from experiments showing that modest uniaxial 
stress is capable of inducing ferroelectricityp 2 ^ Several 
first-principles 22 ' 23 ' 2 ^ 2 ^ 9 , and classical Monte-Carlo sim- 
ulations on an effective Hamiltonian 30 have already pro- 
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vided valuable insight into the structural properties and 
temperature dependence of FE phase transitions in bulk 
STO. In particular, the interaction between FE and AFD 
instabilities in the bulk phase has been thoroughly stud- 
ied from first principles by Sai and Vanderbilt, 31 More- 
over, a previous phenomenological study of the effects of 
epitaxial strain on STO thin films by Pertsev et aZ.p2£ 
based on a Landau theory fit to experimental data from 
bulk phases, produced a rich temperature-strain phase 
diagram and provided support for the idea that the 
oxygen-octahedron rotations have little influence on the 
dielectric response. 

In this work, we compute the effects of the in-plane 
epitaxial lattice-matching constraint on the dielectric re- 
sponse and tunability of SrTi03 using first-principles 
density-functional theory within the local-density ap- 
proximation (LDA). We neglect the AFD instabilities in 
all calculations. We also restrict our analysis to zero tem- 
perature, but for simplicity neglect the quantum fluctu- 
ations. Thus, our focus will be on the ferroelectric soft 
mode and its coupling to strain. Despite these restric- 
tions, the fact that we use a first-principles approach 
means that we do not have to rely on empirical Landau 
parameters as in Ref. l32t Thus, we can confidently make 
predictions under conditions that vary drastically from 
those under which the data determining the parameters 
were obtained, allowing us to consider the effects of large 
epitaxial strains and finite electric fields. 

There is one major limitation of our theory, connected 
with the fact that the LDA tends to underestimate lat- 
tice constants. Because compression tends to suppress 
ferroelectricity, this means that our LDA system is "less 
ferroelectric" than true STO at zero temperature. In- 
deed, we report below a soft-mode frequency of about 75 
cm -1 and a dielectric constant of about 390 at zero tem- 
perature for unstrained STO within the LDA, whereas 
the experimental zero-temperature system is exquisitely 
close to the ferroelectric transition with a; so ft — 10 cm -1 
and e ~ 20, 000. 26 However, raising the temperature also 
has the effect of suppressing the ferroelectric instability, 
so that the experimental system at room temperature is 
characterized by w so f t ~ 90 cm -1 and e ~ 290. Thus, for- 
tuitously, the zero-temperature LDA system corresponds 
reasonably well with the real physical system at a tem- 
perature near, or a bit below, room temperature. Com- 
parisons between these systems must obviously be ap- 
proached with caution, but in fact we find good semi- 
quantitative agreement for several of the physical prop- 
erties of these two corresponding systems as will be pre- 
sented below. 

The manuscript is organized as follows. In Sec. II, we 
provide the details of our approach and approximations 
used, including our handling of finite electric fields. In 
Sec. Ill, we present and discuss the results of our cal- 
culations for epitaxially strained STO, with emphasis on 
the strain dependence of the tunability by finite electric 
fields. Our conclusions appear in Section IV. 



II. METHOD 

A. First principles calculations 

First-principles density-functional calculations are per- 
formed within the local density approximation (LDA) as 
implemented in the PWSCF packaged The exchange- 
correlation energy is evaluated using the Ceperley- Alder 
form with Perdew-Zunger parameterization. Vanderbilt 
ultrasoft pseudopotentialg 35 * are used treating 10 elec- 
trons as valence for Sr (4s 2 4p 6 5s 2 ), 12 electrons for Ti 
(3s 2 3p 6 4s 2 3d 2 ), and 6 for O (2s 2 2p 4 ). To achieve well- 
converged results for small changes in the lattice con- 
stant, the electronic wavefunctions are expanded in plane 
waves up to a kinetic energy of 50 Ry. All calculations 
are performed with a 6x6x6 Monkhorst-Pack k-point 
mesh^i To establish minimum-energy configurations we 
converged the Hellmann-Feynman forces acting on the 
atoms to less than 0.1 mRy/a.u. Density functional per- 
turbation theory (DFPT) is then used to calculate the 
phonon frequencies of the structural optimized systems. 



B. Structural constraints 

The role of epitaxial strain on the structural properties 
of STO is isolated by systematically seeking the ground- 
state structure of five-atom unit cells of bulk STO con- 
strained to several different in-plane lattice constants, 
differing from the theoretical cubic lattice constant by 
fractions ranging between —3% and +3% in steps of 
0.5%, and allowing all atomic positions and the perpen- 
dicular lattice constant to relax fully until the energy 
is minimized. To locate the phase transition points, and 
for the electric field calculations, the in-plane lattice con- 
stant was varied in smaller steps of 0.1% in the paraelec- 
tric region. The presence of the strain necessarily low- 
ers the symmetry of the cubic STO system to tetragonal 
at most; further spontaneous symmetry reduction occurs 
for certain ranges of lattice constant. In the remainder 
of this subsection, we provide specific details of each of 
the structures considered in this work and introduce the 
relevant notation. 



1. Epitaxial strains near zero 

The zero-strain paraelectric phase of SrTiOa has the 
ideal cubic Pm3m perovskite structure, in which the 
octahedral oxygen atoms lie at Wyckoff 3c positions 
(|, 5,0), a single Ti atom lies at the body-centered site 
lb (4, |, |), and the lone Sr cation is at la (0,0,0). Our 
calculations result in a theoretical lattice constant a of 
7.285 a.u., i.e., ~ 1.1 % less than the experimental value 
of 7.365 a.u.; this underestimate is expected when using 
the LDA and is consistent with several previous LDA- 
based studies. 



For small tensile or compressive strains, the structure 
of STO remains centrosymmetric, but the symmetry is 
lowered to tetragonal (space group P4/mmm). In this 
"pseudocubic" phase, the cations sit at Wyckoff positions 
la (0,0,0) for Sr and Id (|, ±, ±) for Ti. There are two 
different Wyckoff positions for oxygen: lc (h, h, 0), which 
will be referred to as Oj_ (in the Ti-0 chain along [001]), 
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i, 0, i) which will be referred to as Oil 



and 2e (0 

and 0\\ y and lie in the same (001) plane as Ti. 



2. Compressive epitaxial strains 
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For large enough compressive epitaxial strains, the 
symmetry of the tetragonal phase is found to be low- 
ered further to non-centrosymmetric PJ^mm. This is a 
ferroelectric tetragonal structure with polarization along 
[001]. The Wyckoff positions are la (0,0, z = 0) for 
Sr, lb (i, |, \ + ATi z ) for Ti, lb (±, i, AOjl,), and 2c 
(^,0, \ + AO|| z ) for oxygens. The presence of a zone- 
boundary instability associated with rotation of the oxy- 
gen octahedra has been discussed for this phase in earlier 
work, 28 i 37 and will not be considered further here. 



3. Tensile epitaxial strains 

Above a critical value of tensile epitaxial strain, STO 
becomes ferroelectric and transforms to an orthorhombic 
Amm2 structure. This structure is non-centrosymmetric, 
exhibiting a nonzero polarization along [110]. The Wyck- 
off positions associated with this phase are 2a (z = 0, 0, 0) 
for Sr, 2b (§ + ATi x , 0, ±) for Ti, and 2b (\ + AO ±x , 0, 0) 
and 4e (j + AO\\ x , j + AOji y , i) for oxygens. The struc- 
ture is illustrated in Fig. ^ The cell-doubling oxygen- 
octahedron rotation expected for this phased will not be 
considered in the remainder of this paper. 



C. Electric field 

1. Field-induced forces 

To describe the dielectric behavior of STO under finite 
dc bias, we must evaluate its properties in the presence 
of a homogeneous electric field. This turns out to be 
subtle; fully first-principles methods for computing the 
behavior of periodic systems in finite fields have only 
recently been developed and are still in their nascent 
stagesi^iiS For STO, we resort to a simple and effective 
approximate technique that consists of first computing 
the Born effective charge tensors Z* for each of the ions 
i, and then adding a term —eZ*£ to their Hcllmann- 
Feynman forces fi, where e is the electronic charge and 
is defined to be positive. The structure is relaxed until 
the total force on each atom is close to zero, i.e., until 
f; = —eZ*£. Similar approximations were discussed by 



FIG. 1: The atomic displacements giving rise to the ferroelec- 
tric Amm2 structure in the tensile-strained region are shown 
in a view of the primitive cubic unit cell projected along [001]. 
On the left: all atoms at their ideal, centrosymmetric posi- 
tions. Arrows indicate displacements of the Ti and oxygen 
atoms. The Sr atoms are fixed at the corners of the cell. 
Right: schematic view of the Amm2 structure (the displace- 
ments are exaggerated for clarity). 

Rabe 40 and then implemented and used by Sai. Rabe and 
Vanderbilt^i this approximation was also recently used 
by Fu and Bellaiche^ 

To be more precise, we ideally would determine the 
structure of STO by minimizing, with respect to first 
wavefunctions ipnk and then to atomic coordinates u^, 
the electric enthalpy T per unit cell, 

T[£) = E KS - flP ■ £, (1) 

where Eks is the internal energy obtained from a 
Kohn-Sham functional (e.g., within the local-density or 
generalized-gradient approximation), 38-42 P = Pi on + P e i 
is the total (ionic plus electronic) macroscopic polariza- 
tion, £ is the electric field, and f2 is the volume of the 
unit cell. Eks and P are explicit functions of atomic 
coordinates and wavefunctions ip n k, with the electric 
field £ entering into T only through multiplication of the 
functional P. (Of course, J 7 , Eks and P can also be 
regarded as implicit functions of £ through the depen- 
dence of the equilibrium values of the "0nk upon £.) By 
a Hellmann-Feynman argument^ the total force acting 
on an ion is 

-^ + n^.£ = f i + Z*e.£, (2) 

where the Z* are the Born effective charges; these also 
depend implicitly on £, even at fixed u^, via the ip n k(£)- 
Here, we neglect the field dependence of and Z*, 
computing these at £ = 0, so that the field enters only 
explicitly as the multiplier of Z*e in Eq. <(5J). This is 
essentially the approximation introduced in Ref. IstI Al- 
though we are able to compute the nonlinear behavior 
of the structural and dielectric properties insofar as they 
arise via field-induced lattice displacements, we neglect 
purely electronic nonlinearities. Thus, the results are rig- 
orously correct only to first order in the fields 
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2. Modeling in a reduced subspace 

In practice, we find that a straightforward minimiza- 
tion of T using the forces of Eq. (J2J) leads to numerical 
instabilities when the system is close to a second-order 
phase transition. This problem arises because the en- 
ergy surface has a very shallow minimum (or competing 
minima and saddle points) and the Hessian matrix be- 
comes poorly conditioned. These numerical difficulties 
can be solved by identifying a subspace that spans, to 
a good approximation, the space of field-induced struc- 
tural distortions, and then parameterizing the energy in 
this subspace. This modeling also helps us better under- 
stand the nonlinear effects of the field on the structure 
and dielectric properties. 

To identify the relevant subspace, we begin by finding 
the pattern of atomic displacements produced by an in- 
finitesimal electric field; it is obtained by multiplying the 
inverse of the force constant matrix i 



>ij,af3 



with the forces 



/i. a = —eJ2p "^i,ap£p- Explicitly, the displacement u,- lC 
of atom i in Cartesian direction a is 



Hi.™ — 



(3) 



In this work we will consider electric fields applied ex- 
clusively along the z direction, perpendicular to the sub- 
strate assuming an (001) film, so the sum over (3 above 
is reduced to a single term. We focus our attention in 
this part of the work on two structures (PA/mmm and 
Amm2); for these (as well as for Pimm, which we will 
not consider further) the symmetry is such that the forces 
and displacements are also only along z. Thus, Cartesian 
indices are dropped for the remainder of this subsection, 
with all quantities referring implicitly to z components 
only. 

Once these displacements Ui have been found, an ex- 
pression for the electric enthalpy T{£) is obtained as fol- 
lows. We gather the Ui into a n-dimensional vector £ 
normalized to unity. (Here n is the number of atoms in 
the unit cell; n = 5 for STO.) For each misfit strain, 
£ defines the subspace of possible field-induced displace- 
ments. The field dependence of the structure and po- 
larization obtained in this model are then that obtained 
by relaxing the ions and minimizing the electric enthalpy 
subject to the constraint that the atomic displacements 
lie along £. This should be a very good approximation 
for small displacements considered here. 

After constructing £, we express the electric enthalpy 
in terms of the scalar amplitude u. Keeping terms only 
to fourth order, we obtain 



F(u) = E a 
= E Q 



bu 2 



bu 2 



du 4 - nAPS 
du 4 — u ^iZ*£. 



structural parameters. For each misfit strain, the ampli- 
tude u was varied up to 0.03 c in steps of 0.002 c (where c 
is the lattice constant perpendicular to the implied sub- 
strate) to obtain the expansion coefficients b and d in 
Eq. dJi£ 

The electric field corresponding to a given u is then 
extracted from the equilibrium condition 



dE 
du 

This leads to 



2bu + 4diT 



£(u) 



2bu + Adu 3 



(5) 



(6) 



Thus, having computed the quantities P, £, and T on 
a mesh of u values, we obtain parametric relations be- 
tween these quantities that can be plotted to reveal the 
dielectric behavior of interest in a numerically stable way. 



D. Dielectric response and tunability 

The dielectric function in the frequency range of the 
optical phonons can be written as the sum of electronic 
and phonon contributions, that is, 



(7) 



(4) 



In most insulating perovskite oxides (ferroelectrics and 
related materials), the electronic contribution is rather 
small (e°° ~ 5) and constant, and the static dielectric 
constant e° = e(0) is typically in the range of 20 — 100, 
so the phonon contribution is expected to dominate. In 
this work, we restrict our focus to the static dielectric 
response and its tunability; we calculate and analyze the 
phonon contribution, neglecting e°°. For the remainder 
of the paper, we drop the superscripts and use the generic 
term dielectric response, even though we only compute 
the phonon contribution. Thus, in what follows, the di- 
electric constant e has the meaning of e ph (0). 

At zero bias, the zero-frequency phonon response is cal- 
culated in a straightforward manner using density func- 
tional perturbation theory (DFPT) to obtain the zone- 
center IR-active phonon modes and their frequencies, and 
using the Berry-phase theory of polarization^ to com- 
pute Born effective charges by finite differences. To eval- 
uate the static response in an applied field, and thus the 
tunability, we use the subspace approach presented in 
Sec. Ill C~2l above, which greatly simplifies our treatment 
while retaining a high degree of accuracy. In the follow- 
ing, both methods are described in detail. 



Zero dc bias 



Here Z* is the mode effective charge given below in 
Eq. @, computed at zero field and for the zero-field 



In the absence of an electric field, the lattice contribu- 
tion to the static dielectric permittivity tensor ep can be 
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written 



jph _ \ - 



47re 2 Z^L^Z, 



c a/3 



^ MM 



^7, 



(8) 



which includes a contribution from each of the zone- 
center polar modes m. Here, f2 is the volume of the 
primitive unit cell, Mq is a reference mass taken as 1 
amu, and Z^ na is the mode effective charge, 




»(*7)> 



(9) 



where e m (i"f) is a dynamical matrix eigenvector. The cor- 
responding real-space eigendisplacement of atom i along 
7 is given by u m (ij) = e m (ij) / \fMl- The Born effective 
charge Z* al {i) is given by 



n dp a 



dui. 



(10) 



In practice, we evaluate the effective charges by finite 
differences, computing the change in polarization AP in- 
duced by several small mode amplitudes Aiti 7 via the 
Berry-phase approach using a 6 x 6 x 20 k-point mesh. 




-1 1 
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FIG. 2: Polarization as a function of misfit strain. Solid 
circles and squares denotes polarization along [001] and [110], 
respectively. 



III. RESULTS AND DISCUSSION 
A. Response to epitaxial strain 

1. Structural properties 



2. Nonzero dc bias 

To calculate the lattice contribution to the static di- 
electric constant in an applied electric field, we use the 
subspace approach presented in Sec. Iff C 21 The field- 
induced change in the structure is specified by a single 
parameter u, which determines the atomic displacements 
through the normalized displacement vector £. We ex- 
press the lattice contribution to the dielectric suscepti- 
bility in terms of the change in polarization induced by 
an applied electric field, 

dP(£) _dP_du 
X ^ ' ~ dE ~~du~d£ 

i 

The expression e = 1 + Airx(£) is used to convert suscep- 
tibility to dielectric constant. In practice, the derivative 
dujdE is calculated numerically once £(u) is determined 
from Eq. © . For strains close to the phase boundary, the 
zero-field dielectric constant computed in this manner is 
identical, by construction, to that obtained in Sec.II.D.l 
using DFPT to obtain phonon frequencies and eigenvec- 
tors, and the Berry-phase calculations to obtain Born 
effective charges. Furthermore, for strains where the pa- 
rameter b was fit to energies for finite u, as described in 
Sec. II. C. 2, the zero-field dielectric constant is in excel- 
lent agreement with the DFPT results. To obtain the 
tunability, we compute the dielectric response for an ap- 
propriate range of values of £(u). 



As described in Sec. Ill Bl we first find the minimum- 
energy structure of SrTiC>3 for values of the misfit strain 
between —3% and +3%. For compressive strains less 
than —0.75%, the lowest-energy structure is ferroelec- 
tric tetragonal PAmm, with polarization along [001]. At 
—0.75%, there is a continuous transition to the nonpolar 
tetragonal PA/mmm phase. At +0.54%, there is another 
continuous transition to the ferroelectric orthorhombic 
Amm2 structure, with polarization along [110]. 

Figure [21 shows the polarization along [001] and [110] 
directions. The polarization increases dramatically with 
strain, and for large strains above 2% (tensile or compres- 
sive), the magnitude of the polarization becomes com- 
parable to that of bulk BaTiC>3, a prototypical ferro- 
electric. This suggests that by simply choosing the ap- 
propriate substrate, the polarization of thin ferroelectric 
STO films, when under short-circuit electrical boundary- 
conditions, could be tuned to a wide range of values. 

In Fig. we show the c/a ratio as a function of misfit 
strain. For tensile (positive) misfit strains, c/a drops 
almost linearly as the magnitude of the misfit strain 
increases and is largely unaffected by the development 
of the in-plane ferroelectric instability at the transition 
to the Amm2 phase. For compressive epitaxial strains, 
in contrast, a noticeable nonlinearity emerges when ap- 
proaching the transition to the PJ^mm phase; this may 
be indicative of the mounting structural frustration that 
is eventually relieved by the occurrence of the transition. 
Once the phase boundary has been crossed, the onset 
and growth of the c-axis polarization in the PJ i mm phase 
further increases the c/a ratio, and c/a once again be- 
comes very nearly linear with misfit strain in the strongly 
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FIG. 3: Calculated c/a as a function of misfit strain in 
SrTiOs. The behaviors for c/a < -1.5% and c/a > +1.5% 
(not shown) are very nearly given by linear extrapolation. 
The open circles denote the samples measured by Hyun and 
Char (Ref.[l3). 



compressive regime. Also shown in Fig. are some mea- 
surements of Hyun and Chasi^ at 77 K that, while noisy, 
show a trend that is roughly consistent with the theory, 
as will be discussed further in Sec. IIII B 21 

It is of interest to compare our misfit phase diagram 
with that of Pertsev et alm££& The latter is obtained by 
expressing the free energy as a function of misfit strain, 
temperature, polarization, and several additional order 
parameters corresponding to the linear oxygen displace- 
ments that account for possible rotations of the oxygen 
octahedra. The parameters in this Landau-theory ex- 
pression are obtained phenomenologically. 

Regarding the nature of the polarization in the se- 
quence of phases at zero temperature, the present first- 
principles results are in very good agreement with the 
Landau analysis. There are several detailed differences, 
however. First, our window of strains over which the sys- 
tem remains paraelectric (between —0.75% and +0.54%) 
is noticeably wider than that previously obtained (be- 
tween -0.2 and -0.02% in Ref. I3I3I . Our overesti- 
mate of the stability of the paraelectric phase can largely 
be attributed to the LDA. Specifically, since our zero- 
temperature calculations within the LDA result in a 
slightly smaller lattice constant than experiment, and 
since the smaller volume stabilizes the paraelectric phase, 
our description of the bulk system yields a more sta- 
ble paraelectric phase, as our underestimate of the bulk 
static dielectric constant attests (calculated ~ 400 com- 
pared to the ~2xl0 4 observed experimentally). In addi- 
tion, octahedral rotations, not included in our analysis, 
might also alter the critical strain corresponding to the 
FE phase transitions in this region. 

An other difference between the present work and 
Refs. I32l33l is the prediction in the latter of a low- 
temperature [100]-polarized ferroelectric phase under in 
a narrow window of tensile strain. Our results are consis- 
tent with only a single phase in this strain regime, with 
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FIG. 4: FE soft-mode frequencies as a function of misfit 
strain. Solid circles indicate the FE soft mode polarized along 
[001]. Solid squares indicate the FE soft mode polarized along 
[110]. The vertical lines indicate the phase transition points 
at —0.75% and +0.54% misfit strain. Open circles denote 
the FE soft-mode frequency if the polarization along [110] is 
suppressed by keeping the atoms on their centrosymmetric 
positions. 



a polarization along [110]. Their polar [100]-oriented or- 
thorhombic phase might be stabilized by the rotations 
not included here, but since it occurs only in a very small 
strain region, we did not investigate it further. For larger 
tensile strains, Pertsev et al£& predicted a phase where 
the polarization is directed along [110] (like ours), with 
an octahedral rotation around the same axis. 



2. Dielectric properties 

In this section, we compute the misfit strain depen- 
dence of the lattice contribution to the dielectric response 
e in zero electric field; finite electric fields will be consid- 
ered in the next section. First, we discuss separately the 
zone-center phonons and the Born effective charges that 
together determine e according to Eq. JHJ. 

The zone-center phonons of each phase are calculated 
at each misfit strain using density-functional perturba- 
tion theory. The lowest-frequency (softest) polar modes 
dominate the dielectric response, and we show the fre- 
quencies of the softest in-plane and out-of-plane trans- 
verse zone-center optical phonons as a function of mis- 
fit strain in Fig. Their behavior reflects the phase 
transition sequence discussed in the previous section. 
By symmetry, these modes are degenerate at zero mis- 
fit strain. However, the lowest-frequency (or soft) mode 
polarized along [001] softens to zero as the misfit strain 
approaches the critical value of —0.75%, signaling the 
second-order transition from the paraelectric PJ^/mmm 
to the [001]-polarized PJ^mm structure. Similarly, the 
lowest-frequency in-plane soft mode softens to zero at 
+0.54% misfit strain, marking the second-order transi- 
tion from P4-/mmm to Amm2, which has its FE polariza- 
tion along the [110] direction. As can be seen in the fig- 
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misfit strain T%1 

FIG. 5: Calculated Born effective charges Z33 (full 
symbols) and Z*i (open symbols) as a function of mis- 
fit strain. The cubic-structure values (at zero misfit 
strain) are Z| 3 (Ti)=7.25, Z 3 * 3 (Sr)=2.55, ZJ 3 (Oi)=-2.06, and 
-£33(03) =— 5.69. For O atoms, up-triangles represent O3, 
down-triangles O2, and diamonds Oi. 



ure, the in-plane soft mode is only weakly affected at the 
P4/mmm-P4mm transition, while the out-of-plane soft 
mode shows a significant hardening in the Amm2 phase. 
The hardening is produced by coupling to the in-plane 
polarization that develops in the Amm2 phase. If the 
polarization is suppressed by keeping the atoms on their 
ccntrosymmetric positions, the mode evolves smoothly 
with increase of the in-plane lattice constant through the 
phase boundary, as is shown by the open circles in Fig.0] 
We now examine the Born effective charge tensors Z* 
as a function of misfit strain. In Fig. we present, for 
brevity, only the Z* x and Z 33 components of the ten- 
sors, where '1' and '3' refer to the [100] and [001] di- 
rections, respectively. While the Born effective charge 
tensors are not diagonal in the Cartesian frame in the 
Amm2 phase (whose principal axes, by symmetry, are 
along [110], [110], and [001]), the computed off-diagonal 
components (not shown) are found to be quite small. In 
what follows, Oi, O2, and O3 refer to the oxygen atoms 
forming Ti-0 chains in the x, y, and z directions, respec- 
tively. 

Figure[S]shows that in the paraelectric region, the Born 
effective charges are very close to the values for the cubic 
structure, 7.25 for Z 3 * 3 (Ti) and -5.69 for Z 3 * 3 (0 3 ). As 
is well known, the fact that these are anomalous (in the 
sense of exceeding the nominal valence) arises from the 
hybridization of Ti and O orbitals in the Ti-0 3 chains 
and is quite sensitive to polar distortions of the chain^ 
Thus, in the ferroelectric phases, we expect a significant 
misfit strain dependence of the Born effective charges for 
these atoms. The effective charge of Ti, Z 33 (Ti), drops 
by almost 16% in the ferroelectric region for compressive 
strain, with a corresponding increase of Z 33 (0 3 ). In com- 
parison, Z*(Sr) is rather insensitive to misfit strain over 
the whole region, showing a weak trend towards lower 
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FIG. 6: Dielectric constant as a function of misfit strain. 
Solid circles and squares denote £33 and en respectively, where 
3 always denotes the [001] direction and 1 refers to [100] in 
the P4mm region and [110] in the Amm2 region. Solid lines 
are fits proportional to (77 — 7/ c ) _1 , where r\ and rj c are the 
actual and critical misfit strains, respectively. Open circles 
are the £33 measured by Hyun and Char (Ref. 



values with increasing misfit strain. 

In the tensile-strain-induced A mm2 ferroelectric phase, 
the orthorhombic symmetry (see Fig. ^) implies that 
Z^Oi) = Z* 2 (0 2 ) ? Z^(0 2 ) = Z* 2 {O x ). In the range 
of strains shown, Zj t 1 (Oi) shows a continuous increase 
of 18%, while the Z\ x for 2 and 3 split slightly. The 
charge neutrality sum rule for the Z* is maintained by a 
corresponding decrease for Z\ x (Ti) . 

In Fig. we show the lattice contribution to the 
static dielectric response in zero electric field over the 
full range of computed misfit strains. The softening to 
zero frequency of the relevant phonons at the second- 
order phase boundaries produces divergences in the di- 
electric response near the critical strains. These transi- 
tions show nearly perfect inverse-power-law behavior, ex- 
cept at the points closest to the phase transitions where 
the low eigenfrequencies lead to numerical inaccuracies. 
Throughout the paraelectric phase, both en and e 33 are 
well over 100. The abrupt drop of e 33 in the tensile-strain 
region is related to the hardening of the lowest-frequency 
phonon polarized along [001] shown in Fig.Q] The com- 
parison with the experimental data of Hyun and Chai^ 
will be discussed further in Sec. IIII B 21 



B. Response to electric field 

Application of a finite electric field to epitaxially 
strained STO leads to changes in the structure and di- 
electric response that depend on misfit strain. In par- 
ticular, the sensitivity to applied field is expected to be 
largest near the phase boundaries discussed above. In 
this section, we focus our attention on the strain regime 
corresponding to the zero-field paraelectric phase, with 
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TABLE I: Normalized displacement vector £ for selected 
strained states in the paraelectric phase. 



strain 




u(0 3 ) 


u(Oi) 


-0.7% 


0.0929 


-0.5081 


-0.6055 


-0.5% 


0.0839 


-0.5111 


-0.6049 


-0.3% 


0.0726 


-0.5129 


-0.6049 


-0.1% 


0.0616 


-0.5143 


-0.6049 


0.0% 


0.0571 


-0.5156 


-0.6045 


0.1% 


0.0518 


-0.5158 


-0.6047 


0.3% 


0.0424 


-0.5176 


-0.6043 


0.5% 


0.0325 


-0.5192 


-0.6039 



particular interest in the behavior as the transition to 
the Pimm phase is approached. This regime is most rel- 
evant to practical applications, as the ferroelectric state, 
its associated hysteresis, and the presence of ferroelectric 
domains are generally undesirable for tunable device ap- 
plications. We also consider electric fields only along the 
[001] direction, in which case all displacements also occur 
only along [001]. 

To briefly re-cap our approach as presented in 
Sec. lII C 2l we parametrize the electric field dependence of 
the atomic displacements up to fourth order. Combined 
with a linearized form of the functional dependence of 
the polarization on atomic displacement, this allows us 
to express the lattice contribution to the dielectric sus- 
ceptibility in terms of the change in atomic displacements 
with electric field. We first discuss the displacements in- 
duced by the applied electric field, and then, in the fol- 
lowing subsection, the resulting lattice contribution to 
the static dielectric response. 



2000 



1500 



| 1000 

a. 

500 



0.0% misfit strain 
-0.1% misfit strain 
-0.2% misfit strain 
-0.3% misfit strain 
-0.4% misfit strain 
-0.5% misfit strain 



-30 



-20 



_j , i i i 

-10 10 
E-Field [Wum] 



20 30 



FIG. 7: Dielectric constant vs electric field for epitaxial STO 
in the paraelectric phase, for a series of compressive strains 
approaching the transition to the ferroelectric PAmm phase. 



very similar to the atomic displacement patterns of the 
[001]-polarized soft mode at zero electric field. For ex- 
ample, the normalized atomic displacements of our com- 
puted normalized soft-mode eigenvector at —0.7% (near 
the critical misfit strain of —0.75%) are (Ti, O3, Oi) = 
(0.0910, -0.5103, -0.6081), almost identical to the dis- 
placement vector £ at the same misfit strain. This shows 
that close to the phase boundary the structural response 
to an electric field is almost entirely dominated by the 
soft mode. 



2. Dielectric properties and tunability 



1. Structural properties 

Using the approximate treatment of electric fields de- 
scribed in Sec. Ill C II we calculate the misfit-strain de- 
pendence of the displacement response to a small elec- 
tric field. The resulting displacements, defined relative 
to the Sr atom, are then normalized to give the vector 
£ reported in Table [I] for selected misfit strains. The O 
displacements in an electric field are relatively insensitive 
to misfit strain, while the Ti displacement in an electric 
field grows with in-plane compressive strain, reflecting 
the fact that in-plane compression, which is accompa- 
nied by a substantial expansion of the c lattice constant 
(see Fig. OH, leads to an opening of the oxygen octahedron 
in the ^-direction and compression in the xy plane. The 
negative sign of the displacements of O atoms in a posi- 
tive electric field results from the negative signs of their 
Born effective charges, while Ti is expected to exhibit 
a positive displacement since its Born effective charge is 
larger than that of Sr. 

For the larger compressive strains in Table HI we find 
that our field-induced displacement vectors (Eq. J3J are 



We now investigate the electric-field dependence of the 
lattice contribution to the dielectric response along the 
[001] direction as the in-plane compressive strain ap- 
proaches the phase boundary with the PAmm phase. The 
results are shown in Fig. [7J As the magnitude of the mis- 
fit strain approaches the critical value of —0.75%, the 
dielectric response at low electric fields grows substan- 
tially. 

We now gain further insight into our calculations 
through a quantitative comparison with experiment. The 
curves in Fig. |7| resemble Lorentzians, and this is ex- 
pected from a phenomenological analysis, as follows. 46 
Within the Landau-Devonshire formalism, the ferroelec- 
tric phase transition can be described by a free-energy 
functional F expanded about the paraelectric phase in 
even powers of the polarization P, i.e., 

F(P, T) = F + AP 2 + BP 4 + CP 6 + ... (12) 

The coefficients depend on misfit strain and temperature; 
A is generally assumed to have the strongest dependence, 
while the variation of B and higher order coefficients is 
small or even negligible. Keeping only terms in F to 
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TABLE II: Zero-field dielectric constant e(0) and field scale 
£o (reported as the slowly varying combination e(0) 3 / 2 £o) ob- 
tained from the fit of Eq. 1161 for each strain state in the 
paraelectric compressively-strained region. 



strain 



e(0) 



<{Qf' 2 £ [V/c 



0.0% 
-0.1% 
-0.2% 
-0.3% 
-0.4% 
-0.5% 
-0.6% 
-0.7% 



391 
473 
584 
760 
1088 
1877 
4125 
19903 



8.3 
9.0 
9.3 
9.6 
10.0 
10.3 
10.1 
10.4 



fourth order, the electric field is then given by 



BP 

£ = — = 2AP + 4BP 3 , 
BP 



and the dielectric susceptibility by 

l_d£_ 

X ~ dP 



= 2A + YIBP 2 



(13) 



(14) 



In the present discussion, the system is in the paraelectric 
phase (A > 0) so that these relations uniquely determine 
functions P{£) and %{£)■ From Eqs. JT3J) and iflijl. it 
is easy to see that in this fourth-order approximation we 
expect x —* constant for P — > 0, while x °c £~ 2 / 3 at large 
P. A useful approximate interpolation formula is then 



x{£) = x(o) 



-1/3 



(15) 



as has been used previously in the literature4£^& Making 
the approximation e 3> 1 so that e ~ 477%, we can write 



e{£) = e(0) 



1 'f 



-1/3 



(16) 



with e(0) = 4ttx(0). 

3/2 

A detailed analysis shows that £q cx e , where the 
constant of proportionality is determined by B, indepen- 
dent of A. That is, e(0) 3 / 2 £o is not expected to depend 
strongly on proximity to the ferroelectric phase transi- 
tion, and thus is expected to vary only slowly with tem- 
perature and misfit strain, and to be only weakly affected 
by the LDA lattice-constant error. Fuchs et aim* used 
Eq. (|16f) to fit their experimentally-measured e-vs.-£ data 
at 200 K, and they observe dielectric constants between 
1480 and 5270; they find values of e(0) 3 / 2 £ ranging from 
7.5 V/cm to 14.3 V/cm for films of 200-500 nm thickness. 
Fitting our theoretical data to the same form, we extract 
values of e(0) 3 / 2 £o ranging from 8.3 V/cm to 10.4 V/cm, 
in general agreement with the findings of Fuchs et al. The 
results of our comparison are summarized in Table ITU 



TABLE III: Measured misfit strain, c/a, and zero electric 
field dielectric constant for STO films on various substrates. 14 
CRO and LAO stand for CaRuOa and LaAlOa, respectively. 



strain 


c/a 


e(£ = 0) 


materials 


-0.3073% 


1.0072 


640 


Au/STO/CRO/LAO 


-0.1536% 


1.0056 


690 


Au/STO/CRO/STO 


+0.2305% 


0.9937 


400 


Au/STO/SRO/LAO 


+0.7170% 


0.9942 


360 


Au/STO/SRO/STO 



Hyun and Char— have grown SrTiC>3 on different sub- 
strates to investigate the influence of epitaxial strain on 
the tunability of SrTiOa. For four samples, they report 
measurements at 77 K of in-plane lattice constants, c/a 
ratios, and dielectric constants, as summarized in Table 
Till There is good agreement between measured and cal- 
culated c/a ratios for the values of misfit strain observed 
in Ref. IliL as shown in Fig. 3. For three of the sam- 
ples, the dielectric constant increases with compressive 
misfit strain in the paraelectric phase, in agreement with 
our calculations. The first sample is an exception: we 
expect it to exhibit a higher tunability, given its strain 
state. This discrepancy can be attributed to the poorer 
quality of samples grown on LaAlOa^* Comparing our 
zero-temperature calculations in IIII(I with the measured 
value of 690 for the —0.1536% sample, and an additional 
computation of 319 for strain = +0.2% with the mea- 
sured value of 400 for the +0.2305% sample, (see Ta- 
ble ^HJ, we find surprisingly good agreement. This may 
be partly due to slight shifts in the experimental phase 
boundaries owing to finite temperature; because of our 
underestimate of lattice constants, our paraelectric phase 
is over-stabilized, mimicking the effect of temperature 
and shifting the calculated phase boundary towards more 
compressive misfit strains. 

Unlike the first three samples listed in Table ITTTI which 
we expect to be paraelectric based on their in-plane lat- 
tice constants, the fourth sample should already be in the 
orthorhombic, tensile FE region, according to our calcu- 
lations and to the phase diagram of Pertsev et al. For 
this reason, we did not include this experimental data 
point in Fig. |SJ We should like to note, however, that 
the trends in the dependence of c/a and e(0) on in-plane 
strain in Table ITTTI do not seem consistent regarding this 
last data point. As the strain changes from +0.23% to 
+0.72%, the measured c/a and e(0) remain almost iden- 
tical. The behavior of c/a, for example, seems inconsis- 
tent with the roughly linear dependence of c/a on misfit 
strain predicted in our Fig. EI suggesting that more thor- 
ough checks of the experimental behavior in the region 
of strong tensile strain may be called for. 

Finally, we have calculated the strain dependence of 
the tunability parameter in the paraelectric phase for 
two different values of the bias field, £ = 10 V/ /im and 
£ — 25 V//im, with results shown in Fig. [21 Although 
a direct comparison is not possible because the theoreti- 
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FIG. 8: Tunability n = [e(0) - e(£)]/e(0) as a function of 
misfit strain. The values for e(£) were taken at £ = 10 V//xm 
(full squares), and £ — 25 V/fim (full circles). Open circles 
are the experimental results of Hyun and Char for £ = 10 
V//xm. 

cal and experimental conditions are rather different, the 
measurements of Hyun and Char for £ = 10 V//im are 
also shown. These results are closer to our calculated 
tunability for £ = 25 V/^m than the one for £ = 10 
V//xm, consistent with the fact that our LDA-computed 
zero-temperature system is further from its ferroelectric 
transition (as indicated, e.g., by the smaller zero-field di- 
electric constant) than for the experimental system at 
77 K, and therefore less tunable. Nevertheless, our re- 
sults do tend to confirm that a substantial variation in 
the tunability can be attained by growing STO on sub- 
strates having different lattice constants. 



first-principles density-functional theory within the local 
density approximation. We computed the tunability of 
STO, and found this quantity to be highly sensitive to 
epitaxial strain. We also find that the dielectric constant 
itself varies significantly with strain, and for sufficiently 
large compressive strains, STO can be made ferroelectric 
with a polarization similar to that of bulk BaTi03. Our 
results are in good agreement with available experiments. 

Our studies complement and extend previous theo- 
retical work on STO films using a phenomenological 
Landau theory. Our parameter-free description pro- 
vides structural parameters and phonon frequencies that 
can be compared with future experiments; in addi- 
tion, we observe significant differences between our zero- 
temperature phase diagram and that obtained using Lan- 
dau theories fit to empirical data. 

These results should be useful for the analysis of ul- 
trathin epitaxial films and superlattice structures in- 
volving STO. In the latter, high strain states, achieved 
through layering, can result in novel artificial materials 
with enhanced properties over those of their bulk con- 
stituents, as was recently predicted for BaTiOa/SrTiOs 
superlattices4^ 

For coherent epitaxial systems more complex than 
STO and at finite temperatures, it seems clear that epi- 
taxial strain will generally be a crucial factor in deter- 
mining properties. In particular, we expect our results 
for STO to be a useful guide in understanding the strain 
dependence of more complex ferroelectric thin films with 
high tunability, such as BST. 



IV. SUMMARY 
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